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A method of monitoring a magnetic field drift of a magnetic resonance imaging apparatus 



FIELD OF THE INVENTION 

The present invention is related to the field of magnetic resonance imaging 
(MRI), and more particularly without limitation to magnetic field drift compensation. 

5 BACKGROUND AND PRIOR ART 

Prior art magnetic resonance imaging apparatus' generate a uniform magnetic 
field which is also referred to as the main field or the polarising field B 0 . The purpose of the 
magnetic field is that the individual magnetic moments of the spins in the tissue to be 
visualised attempt to align with this polarising field, but precess about it in random order at a 
10 characteristic Larmor frequency which is determined by the gyromagnetic constant of the 
spins and the polarising magnetic field Bo. 

However during operation of a magnetic resonance apparatus the magnetic 
field can vary over time due to various reasons. For example when the magnetic resonance 
apparatus is under heavy duty the passive iron shim pieces of the coils which generate the 
IS magnetic field heat up which leads to a fluctuation of the magnetic field. 

U.S. Patent No. 6,294,913b 1 shows a method for compensation of variations 
in the polarising magnetic field during magnetic resonance imaging. Monitor signals are 
acquired in an interleaved manner during a scan with the MRI system. Frequency changes 
caused by variations in the polarising magnetic field B 0 are measured using the monitor 
20 signals, and these measured frequency changes are employed to compensate image data 

acquired during the scan. This compensation is achieved by changing the frequency of the RF 
transmitter and receiver to offset the effects of changes in Bo. A disadvantage of this 
compensation method is that special monitor signals are required for the compensation. 

It is therefore an object of the present invention to provide for an improved 
25 method of monitoring and compensating a magnetic field drift 



SUMMARY OF THE INVENTION 

The present invention provides for a method of monitoring a magnetic field 
drift of a magnetic resonance imaging apparatus which does not require special monitor 
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signals or monitor sensors. Rather the method of the invention can be performed on the basis 
of the magnetic resonance signals which are acquired during magnetic resonance image data 
acquisition, such as by single shot EPI or by a gradient echo sequence. 

In essence the phases of at least two magnetic resonance signals are acquired 

5 an echo time after the corresponding RF excitations. This corresponds to the central k-space 
line which has frequency encoding but no phase encoding. This central k-space is usually 
designated as k y = 0. The difference of two consecutive phase measurements which are 
acquired at a certain time interval provides the shift of the resonance frequency. This enables 
monitoring of the shift of the resonance frequency and compensation of the magnetic field 

10 drift. 

In accordance with a further preferred embodiment of the invention a single 
shot EPI method is used for the magnetic resonance data acquisition. Single shot EPI has the 
advantage that every data acquisition contains the k y = 0 line such that for each data 
acquisition the required phase information can be obtained. 

IS In accordance with a further preferred embodiment of the invention a gradient 

echo sequence is used for the magnetic resonance image data acquisition. During such a 
gradient echo sequence the k-space is scanned along an arbitrary trajectory which typically 
also contains the k-space line k y = 0. For example a full gradient echo sequence has 256 data 
acquisitions one of which is representative of k y = 0. The data acquisition for k y = 0 enables 

20 to obtain the required phase information once for each complete gradient echo sequence. 

In accordance with a further preferred embodiment of the invention extra 
scans along k y = 0 are performed during a gradient echo sequence in order to obtain multiple 
phases during the complete sequence. This enables to monitor a magnetic field drift which 
occurs during a single gradient echo sequence. 

25 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the following preferred embodiments of the invention will be described in 
greater detail by making reference to the drawings in which: 

Figure 1 shows a block diagram of a magnetic resonance apparatus with 
30 magnetic field drift compensation, 

Figure 2 is illustrative of the phase determination for a single shot EPI, 

Figure 3 is illustrative of a flow chart for monitoring and compensating the 
magnetic field drift in the case of single shot EPI, and 
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Figure 4 is illustrative of a flow chart for magnetic field monitoring and 
compensation if a gradient echo sequence is employed 

DETAILED DESCRIPTION 
5 Figure 1 shows a magnetic resonance device 1 which includes a first magnet 

system 2 for generating a steady magnetic field, and also several gradient coils 3 for 
generating additional magnetic fields having a gradient in the X, Y, Z directions. The Z 
direction of the co-ordinate system shown corresponds to the direction of the steady magnetic 
field in the magnet system 2 by convention. The measuring co-ordinate system x, y, z to be 

10 used can be chosen independently of the X, Y, Z system shown in FIG. 3. The gradient coils 
are fed by a power supply unit 4. An RF transmitter coil 5 serves to generate RF magnetic 
excitation pulses and is connected to an RF transmitter and modulator 6. 

A receiver coil is used to receive the magnetic resonance signal generated by 
the RF field in the object 7 to be examined, for example a human or animal body. This coil 

15 may be the same coil as the RF transmitter coil 5. Furthermore, the magnet system 2 encloses 
an examination space which is large enough to accommodate a part of the body 7 to be 
examined. The RF coil 5 is arranged around or on the part of the body 7 to be examined in 
this examination space. The RF transmitter coil 5 is connected to a signal amplifier and 
demodulation unit 10 via a transmission/reception circuit 9. 

20 The control unit 1 1 controls the RF transmitter and modulator 6 and the power 

supply unit 4 so as to generate special pulse sequences which contain RF pulses and 
gradients. The phase and amplitude obtained from the demodulation unit 10 are applied to a 
processing unit 12. The processing unit 12 processes the presented signal values (also 
referred to as k-space) so as to form an image by transformation. This image can be 

25 visualized, for example by means of a monitor 13. 

Processing unit 12 determines a shift A/" of the resonance frequency on the 
basis of the phase information obtained from the demodulation unit and provides the shift 6f 
to the control unit 11 for compensation of the shift 

The acquisition of the required phase information is explained in more detail 

30 in the following by making reference to figure 2. First an RF magnetic excitation pulse 200 
is generated by the transmitter coil (cf. transmitter coil 5 of figure 1). Further a slice selecting 
gradient G z pulse 202, a bipolar switched gradient 204, and a phase encoding gradient 206 
are applied The application of these gradients 202, 204 and 206 for single shot EPI magnetic 
resonance data acquisition is as such known from the prior art. 
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The single shot EPI provides measurement signal 208. The peak of 
measurement signal 208 is the echo time TE after the RF magnetic excitation pulse 200. The 
phase of measurement signal 208 at time TE with respect to the RF magnetic excitation pulse 
200 is determined for the purpose of monitoring a drift of the magnetic field and a 
5 corresponding drift of the resonance frequency. In order to determine the resonance 
frequency drift, if any, the phase of the measurement signal 208 is determined again at a 
subsequent single shot EPI. The phase difference between the single shot EPIs enables to 
calculate the shift of the resonance frequency. This will be explained in greater detail by 
making reference to figure 3. 
10 The phase information at time TE can be obtained from the measurement 

signal 208 in the time domain or in the frequency domain which requires to perform a Fourier 
transformation on measurement signal 208. The required phase information corresponds to a 
central line of the k-space which is usually referred to as k y = 0, i.e. a scan with frequency 
encoding but without phase encoding. The fact that there is no phase encoding enables to use 
IS drifts of the phase information obtained from the measurement signal 208 for determination 
of the resonance frequency shift. 

Figure 3 shows a flow chart for monitoring and compensating of a magnetic 
field drift in the case of single shot EPI. In step 300 the index d is set to zero. The shift of the 
resonance frequency Af is also set to zero in step 300. 
20 In step 302 the first single shot EPI is performed and the measurement of the 

phase <}> a (d) is performed as explained with respect to figure 2. 

The time interval of TE after the first single shot EPI is performed in step 302 
the second single shot EPI is performed in step 304. Again the phase <J) a is measured 

In step 306 a shift Af of the resonance frequency is calculated by calculating 
25 the difference between the phases (j> a which have been acquired in steps 302 and 304 and by 
dividing the difference by 2n x the time interval between the consecutive single shot EPIs. In 
this case this time interval is the echo time TE. 

In step 308 the absolute value of the shift 4f is compared to a threshold value. 
If the threshold value is surpassed a drift compensation is performed in step 310. Preferably 
30 this is done by adjusting the frequency of the RF magnetic excitation pulses to the shifted 
resonance frequency. 
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If the shift 6f does not surpass the threshold no drift compensation is 
required. In step 312 the index d is incremented and the control goes back to step 304 for 
continuous monitoring the shift Af . 

Figure 4 shows an alternative embodiment which uses a gradient echo 
5 sequence rather than single shot EPI. Step 400 corresponds to step 300 of figure 3. In step 
402 the magnetic resonance data acquisition is started by starting to scan the k-space along a 
given trajectory. The trajectory may or may not cover the central line of the k space which is 
ky = 0. 

In step 404 the scan of the k-space along the trajectory is interrupted in order 
10 to perform a data acquisition for the central k-space line with frequency encoding but no 
phase encoding, i.e. k y = 0. This way the phase information <j> a is obtained. 

In step 406 the scan along the k-space trajectory is continued. After a certain 
time interval the scan along the k-space trajectory is interrupted again in order to perform 
another data acquisition for ky = 0 to obtain another phase information <j> a . 
IS In step 410 the shift of the resonance frequency is calculated by calculating the 

difference between the phases determined in steps 404 and 408 and by dividing the difference 
by2rcTE. 

In step 412 the absolute value of the frequency shift Af is compared to a 
threshold value. If the threshold value is surpassed drift compensation is performed in step 

20 414. In step 416 the index d is incremented and the control goes back to step 406 in order to 
continue the scan along the k-space trajectory. This procedure continues during the entire 
gradient echo sequence for continuously monitoring the shift of the resonance frequency and 
compensating the drift of the magnetic field if necessary. 

The RF pulses of this modified gradient echo MRI sequence are applied at a 

25 single frequency (for a single slice); the phase measured at the echo (without encoding 

gradients) reflects only the phase difference accumulated between the RF pulse and the TE in 
case that the NMR resonant frequency (due to B0) is not equal to the excitation frequency 
(excluding susceptibility and chemical shift effects). When the magnet B0 is such that the 
NMR resonant frequency is the same as the excitation frequency of the RF pulse then the 

30 phase error accumulated between RF pulse and TE will be zero or constant from TR to TR. 

When the B0 is such that it corresponds to a different resonant frequency than 
the one the RF pulse is exciting, then the phase error accumulated between RF pulse and TE 
will be a finite value (proportional to the difference and TE). At the peak of each RF pulse, 
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all magnetization is in phase, irrespective of the BO. Only after this time does the BO have an 
effect on die magnetization. 

So, each RF pulse behaves like a phase reset (as fer as the fO measurement 
process is concerned). The relative phase error (relative to zero at the RF pulse peak) at the 
5 TE increases as BO moves further from the RF excitation frequency. Other mechanisms 
(chemical shift, susceptibility) can contribute to the phase error measured at the TE. 
Assuming that these other contributors to phase error are constant from TR to TR, the change 
in BO can be determined by calculating the difference in phase errors between data acquired 
form two different TR's. 

10 
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CLAIMS: 



1 . A method of monitoring a magnetic field drift of a magnetic resonance 

imaging apparatus, the method comprising the steps of. 

- performing a first data acquisition by a first magnetic resonance signal being caused 
by a first excitation, 

5 - determining a first phase of the first magnetic resonance signal an echo time after the 
first excitation, 

- performing a second data acquisition by a second magnetic resonance signal a time 
interval after the first data acquisition, the second magnetic resonance signal being 
caused by a second excitation, 

10 - determining a second phase of the second magnetic resonance signal the echo time 
after the second excitation, 

- determining a shift of a resonance frequency based on a difference of the first and 
second phases. 

IS 2. The method of claim 1, whereby the first and second data acquisition are 

performed using a signal shot EPI method. 

3. The method of claim 1, whereby the first and second data acquisitions are 
performed by means of a gradient echo sequence method. 

20 

4. The method of claim 3, whereby a k-space is scanned and second data 
acquisitions are performed intermittently to determine second phases in order to continuously 
monitor the shift of the resonance frequency. 

25 5. The method of claim 4, whereby the second data acquisitions are performed 

after fixed time intervals. 
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6. The method of any one of the preceding claims 1 to 5, further comprising 

compensating the magnetic field drift by changing the frequency of the excitation in 
accordance with the shift of the resonance frequency. 

5 7. The method of any one of the preceding claims 1 to 6, further comprising 

compensating the magnetic field drift by adjusting the magnetic field. 

8. The method of any one of the preceding claims 1 to 7, further comprising 
comparing the shift of the resonance frequency to a threshold value and compensating the 

10 magnetic field drift if the threshold value is surpassed. 

9. The method of any one of the preceding claims 1 to 8, whereby the first and 
second phases are determined in the time domain. 

IS 10. The method of any one of the preceding claims 1 to 8, further comprising 

performing a Fourier transformation of the first and second magnetic resonance signals and 
determining the first and second phases in the frequency domain. 

11. A computer program product, in particular digital storage medium, for 

20 monitoring a magnetic field drift of a magnetic resonance imaging apparatus, the computer 
program product comprising program means being adapted to perform the steps of: 

- determining a first phase of a first magnetic resonance signal, an echo time after a 
first excitation, 

- determining a second phase of a second magnetic resonance signal the echo time after 
25 a second excitation, whereby the second magnetic resonance signal is acquired a time 

interval after the first magnetic resonance signal, 

- calculating a shift of a resonance frequency based on a difference of the first and 
second phases. 



30 12. The computer program product of claim 1 1, the program means being adapted 

to continuously monitor the shift of the resonance frequency. 
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13. The computer program product of claim 11 or 12, the program means being 

adapted to control an excitation synthesiser in accordance with the shift of the resonance 
frequency. 

5 14. The computer program product of claim 1 1, 12 or 13, the program means 

being adapted to control the magnetic field in accordance with the shift of the resonance 
frequency. 

15. A magnetic resonance imaging apparatus comprising processing means (12) 
10 for determining a first phase of a first magnetic resonance signal an echo time (TE) after a 

first excitation, for determining a second phase of a second magnetic resonance signal the 
echo time after a second excitation, the second magnetic resonance signal being acquired a 
time interval after die first magnetic resonance signal, and for calculating a shift of a 
resonance frequency based on a difference of the first and second phases and the time 
IS interval. 

16. The magnetic resonance imaging apparatus of claim IS having display means 
(13) for displaying of the shift of the resonance frequency. 

20 17. The magnetic resonance imaging apparatus of claims IS or 16 further 

comprising control means (1 1) for controlling the generation of the excitations in accordance 
with the shift of the resonance frequency. 

18. The magnetic resonance imaging apparatus of claims 15, 16 or 17, further 

25 comprising control means for controlling of the magnetic field in accordance with the shift of 
the resonance frequency. 
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ABSTRACT: 



The invention enables to monitor a magnetic field drift of a magnetic 
resonance imaging apparatus on the basis of the magnetic resonance signals, which are 
acquired during magnetic resonance image data acquisition, such as by single shot EPI or by 
a gradient echo sequence. The phases of at least two magnetic resonance signals are acquired 
5 an echo time after the corresponding RF excitations. This corresponds to the central k-space 
line, which has frequency encoding but no phase encoding. The difference of two 
consecutive phase measurements, which are acquired at a certain time interval provides the 
shift of the resonance frequency. This enables monitoring of the shift of the resonance 
frequency and compensation of the magnetic field drift. 

10 

(Figure 3) 
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